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ABSTRACT 



Of the over 700 exoplanets detected to date, over half are on non-circular 
orbits, with eccentricities as high as 0.93. Such orbits lead to time-variable 
stellar heating, which has implications for the planet's atmospheric dynamical 
regime. However, little is known about this dynamical regime, and how it may 
influence observations. Therefore, we present a systematic study of hot Jupiters 
on highly eccentric orbits using the SPARC/MITgcm, a model which couples a 
three-dimensional general circulation model with a plane-parallel, two-stream, 
non-grey radiative transfer model. In our study, we vary the eccentricity and 
orbit-average stellar flux over a wide range. We demonstrate that the eccentric 
hot Jupiter regime is qualitatively similar to that of planets on circular orbits; 
the planets possess a superrotating equatorial jet and exhibit large day-night 
temperature variations. We show that these day-night heating variations induce 
momentum fluxes equatorward to maintain the superrotating jet throughout its 
orbit. As the eccentricity and/or stellar flux is increased, the superrotating jet 
strengthens and narrows, due to a smaller Rossby deformation radius. For a select 
number of model integrations, we generate full-orbit lightcurves and find that the 
timing of transit and secondary eclipse viewed from Earth with respect to periapse 
and apoapse can greatly affect what we see in infrared (IR) lightcurves; the peak 
in IR fiux can lead or lag secondary eclipse depending on the geometry. For 
those planets that have large day-night temperature variations and rapid rotation 
rates, we find that the lightcurves exhibit "ringing" as the planet's hottest region 
rotates in and out of view from Earth. These results can be used to explain 
future observations of eccentric transiting exoplanets. 



-3- 
1. Introduction 

Since the first confirmations in the mid-1990s (Mayor & Queloz 1995; Wolszczan 1994) 
the detection and characterization of extrasolar planets continue to be major fields in 
astronomy and planetary science. Over 700 planets have been detected from the ground 
and space, more than half of which are classified as "hot Jupiters" , Jovian-mass planets 
that orbit their parent stars at distances less than 0.1 AU. A number of these hot Jupiters 
transit their host star along our line of sight, allowing us to observe them as they pass in 
front and behind their parent star (e.g. Charbonneau et al. 2008; Swain et al. 2008; Pont 
et al. 2008; Knutson et al. 2007, 2008, 2009). Using these observations, we can infer much 
about these planets' atmospheric composition, temperature structure, and circulation. 

A fifth of these transiting exoplanets have eccentricities greater than 0.1, with values 
as large as 0.93 (HD80606b, Naef et al. 2001). These eccentric planets are subject to 
highly time- variable heating which has a large effect on the planet's atmospheric dynamics. 
Among planets amenable to observational follow-up, HAT-P-2b, which has an eccentricity 
of 0.52, undergoes a factor of 9 variation in flux throughout its orbit. HD 17156b (e = 0.67) 
experiences a factor of 27 variation in stellar flux, and HD 80606b, with its large eccentricity 
(e = 0.93), undergoes an impressive factor of 828 variation in flux. Because they are 
transiting, we can probe their atmospheres as we can for planets on circular orbits. However, 
while it is clear the strongly variable heating leads to a vastly different forcing regime than 
for exoplanets on circular orbits, it remains unknown whether this causes a fundamentally 
different dynamical regime: is the circulation quantitatively similar to that on circular hot 
Jupiters, or is it a completely new circulation regime? 

Eccentric transiting exoplanets present unique challenges when one attempts to extract 
information about their atmospheres from observational data. In particular, interpretation 
of flux maxima and minima in infrared lightcurves can be complicated by the convolution of 
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spatial effects (for example, hot spots on the planet that rotate into and out of view along 
our line of sight) with temporal effects (planet getting colder /warmer at apoapse/periapse 
passage). Langton and Laughhn (2008) and Cowan and Agol (2011) address this problem as 
applied to particular targets, but use only a two-dimensional and one-dimensional dynamical 
models, respectively. To fully capture these effects, a three-dimensional circulation model 
that self consistently calculates heating and wind velocities is needed. 

Hence, it is crucial to conduct a comprehensive study that establishes the dynamical 
regime, temperature structure, and observational implications of eccentric exoplanets. We 
use a three-dimensional atmospheric circulation model coupled to a non-grey radiative 
transfer scheme to study eccentric hot Jupiters as a whole. In Section 2 we will describe 
our model setup and integrations. Section 3 describes the dynamical regime and its 
dependence on eccentricity and mean stellar flux. Section 4 presents synthetic lightcurves 
and attempts to determine what can be learned about the atmospheric circulations from 
remote measurements. Finally, Section 5 concludes and compares results to known eccentric 
hot Jupiters. 

2. Model 

We adopt the Substellar and Planetary Atmospheric Radiation and Circulation 
(SPARC) model, which couples the MITgcm (Adcroft et al. 2004) with a two-stream 
implementation of the multi-stream, non-grey radiative transfer scheme developed by 
Marley & McKay (1999). A brief discussion of the model is included here; for a more 
complete description, see Showman et al. (2009). 

The MITgcm solves the primitive equations, a simplification of the Navier-Stokes 
equations where the horizontal flow length scale exceeds vertical length scales; these 
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equations are valid in stably stratified atmospheres. Hot Jupiters generally satisfy these 
criteria, with horizontal and vertical length scales of ^ 10^ km and ^300 km, respectively, 
and highly stratified atmospheres. The radiative transfer code solves the two-stream 
radiative transfer equations, and employs the correlated-fc method (e.g. Mlawer et al. 1997; 
Goody et al. 1989; Fu & Liou 1992; Marley & McKay 1999) to solve for upward/downward 
fluxes and heating/cooling rates through an inhomogeneous atmosphere. This method 
retains most of the accuracy of full line-by-line calculations, while drastically increasing 
computational efficiency. The correlated-fc method maps the absorption coefficients, 
denoted by k{iy) where u is frequency, from spectral space {k = A;(z/)) to a space where 
k varies with a variable g {k — k{g)). The inversion of this function, g{k)^ is commonly 
referred to as the cumulative distribution. 



where fi{k) is the distribution function for the absorption coefficient in the zth atmosphere 
layer in a spectral interval Az/. In this space, g varies monotonically with k and k{g) 
can be easily approximated. The correlated-A; method splits g into several subintervals 
each containing 10^ — 10^ values of k[g)^ and determining an equivalent monochromatic 
absorption coefficient, Ki^ for each subinterval. Each value of Ki is used to calculate the 
transmission, T, as a function of column mass, by the equation 



where Wi are the weights corresponding to width of each subinterval. In this way, multiple 
opacity sources can be modeled, accounting for optical effects such as clouds, aerosols, 
dust, scattering and pressure-induced absorption. The fluxes for each of these effects are 
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(2) 
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obtained, summed, and weighted to obtain the total flux in each spectral interval for each 
atmospheric layer, and then further summed to obtain the total, wavelength-integrated net 
upward/downward flux. We then calculate the heating rate by flnite-differencing the fluxes 
and pressures between interfaces over and underlying a given dynamical level, 

dF 

« = Pa^ (3) 

Here q is the heating rate, g is the planet's gravity, and ^ is the gradient of flux with 
pressure. For each simulation, we utilize a cubed-sphere grid with a horizontal resolution 
of C32 (approximately equivalent to 64x128 in latitude and longitude) and Nl =40 or 76 
pressure levels. The lowermost Nl — 1 levels extend from a mean pressure of 200 bars at 
the bottom to 0.2 mbar at the top, evenly spaced in log pressure. The top level extends 
from a pressure of 0.2 mbar to zero. 



2.1. Updates to SPARC/MITgcm 

2.1.1. Reduced number of frequency bins 

While the models of HD 189733b and HD 209458b in Showman et al. (2009) used 
30 frequency bins to model circulation and heating on each planet, here the opacities are 
statistically weighted into 11 frequency bins to improve computational efliciency. For both 
the 30- and 11-bin implementations, the radiative transfer is calculated at 8 values of k{g) 
in each of the frequency bins. Four of the values of k{g) are sampled at < 0.95, meaning 
that four calculations are done for the weakest 95% of the spectral lines. The other four 
radiative transfer calculations are done for the strongest 5% of the lines {g > 0.95). This 
allows us to calculate the radiative transfer for both strong and weak spectral lines. 
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We list the 11 new frequency bins in Table 1. Before implementing the updated 
scheme into our simulations, we conducted one-dimensional and three-dimensional tests to 
reproduce the results for HD1897833b from Showman et al. (2009). Once successful, we 
applied the new scheme to the model integrations described here. (Appendix A describes 
these validation tests in detail.) 

2.1.2. Updated collision-induced absorption due to H2 — H2 and H2 — He collisions 

The coUision-induced absorption (CIA) by H2 molecules due to H2 — H2 and 
H2 — He collisions has recently been recomputed from first principles (Saumon et al. 
2012 and references therein). The new formulation uses state-of-the-art quantum 
mechanical calculations of the potential energy surface and of the interaction-induced 
dipole moment surface for colliding H2 — H2 and H2 — He pairs. This new calculation 
is more reliable at higher temperature and higher photon energies than the earlier 
calculations by Borysow and collaborators that radiative transfer groups have been using 
(http://www.astro.ku.dk\^aborysow/) . We have included these new tabulations of H2 
CIA opacity in our code. 

2.2. Grid simulations 

Table 2 lists our model integrations for this study. For all models, we assume the 
radius and gravity of a generic hot Jupiter^. Figure 1 shows the parameter space plotted 
as a function of semi-major axis and eccentricity, along with detected exoplanets with 

^HD189733b, Rp = 8.2396 x 10^ m, g = 21.4 m s-^ Note that we could have easily 
chosen any other well-characterized hot Jupiter, such as HD 209458b. However, the primary 
goal of this study is to illustrate the dynamical and observational effects of eccentricity and 
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semi-major axes less than 0.5 AU. In these simulations, we systematically vary both the 
average stellar flux and eccentricity. Orbital eccentricities of 0.0 to 0.75 and orbit-averaged 
stellar fluxes, (F) of 468183, 185691, 73680 and 11617 Wm-^ were explored. These values 
correspond to equilibrium temperatures (Teq) of approximately 1199, 951, 755 and 476 
K, respectively. For these given values of e and (F), the orbital semi-major axis, a, is 
calculated using the equation 



where L is the luminosity of HD 189733. These chosen simulations fall well within the 
current exoplanet population. 

We assume the atmospheric composition to be 1 x solar metallicity without TiO and 
VO; this composition has had success in explaining observations of HD 189733b. As in 
Showman et al. (2009), opacities are determined assuming local chemical equilibrium 
(accounting for rainout of condensates) at each temperature and pressure over the 3D 
grid, assuming the elemental abundances of Lodders (2003). We set Cp and k to values 
appropriate for a i/2-dominated atmosphere (1.3 x 10^ J kg~^ and 2/7, respectively). 

For eccentric cases, we assume the planet is pseudo-synchronously rotating its host 
star — we assume the planet's tidal interactions with the star force the same side of the 
planet to approximately face the star every periapse passage. We calculate this using the 
Hut (1981) formulation: 



L 



(4) 



47ra2(l -e2)i/2- 




(5) 



rot -^or6 



2 ^ ^8 ^ 16 



average stellar flux. The precise planetary radius and gravity are of secondary importance. 
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In each integration, we assume the winds to be initially zero, with each vertical 
atmospheric column assigned the same global-mean, radiative equilibrium pressure- 
temperature profile calculated using the same opacities as the main model (Figure 2). 
These temperature-pressure profiles are calculated for each value of average stellar flux (see 
above) using the methods described in Fortney et al. (2005, 2008). 

Most model integrations were run until the winds reached a steady state. Figure 3 
shows the root-mean-squared (RMS) velocity plotted as a function of pressure and simulated 
time, for a model integrated for >6000 Earth days. While the winds are not fully in a 
steady state, drag could be artificially implemented to force steady state winds. 

3. Results 

3.1. Circulation regime 

To illustrate the circulation regime of an eccentric exoplanet, we present a nominal 
model where e = 0.50, a = 0.0534 AU, and T^q = 951 K, representative of most of the 
simulations (Figure 4). The top panel shows the horizontally averaged global temperature 
as a function of pressure, time (in Earth days) and true anomaly, /, where / = 0° at 
periapse and / = ±180° at apoapse. Temperatures remain fairly constant at pressures 
exceeding 1 bar. However, in the uppermost layers of the atmosphere, temperatures vary 
throughout the orbit, peaking ^6 hours after periapse passage. The bottom panel shows 
the zonal-mean zonal wind^ of the same model time- averaged over a full orbit. The flow 
exhibits a superrotating (eastward) equatorial jet which is the dominant feature of the flow, 
with peak winds of about 5000 ms~^. 

^ Zonal wind is defined as east- west winds with positive (negative) values denoting east- 
ward (westward) winds; a zonal average is an average in longitude. 
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Snapshots of the wind and temperature profile throughout the orbit at a pressure level 
of 30 mbar (approximately the level of the infrared photosphere) show that the planet 
undergoes a flash heating event as the planet approaches periapse passage (Figure 5). The 
eight panels show snapshots throughout one complete orbit from apoapse (top figure) to 
hours after periapse (/ = 119°, follow arrows) and back again to apoapse, all with the same 
colorscale for temperature. Like the circulation of HD 189733b in its nominally circular 
orbit (see Showman et al. 2009), the flow maintains the eastward equatorial superrotating 
jet throughout its orbit (Figure 6). As the planet approaches periapse (Figure 5, / = —12°), 
temperatures rise to ^1000 K at the equator, and a high-amplitude eddy structure develops 
with a pattern of temperature and wind vectors that tilt eastward toward the equator. 
This so-called "chevron" -shaped feature is prominent at periapse (bottom panel), extends 
over a wide range of latitude and longitude, and persists hours after periapse. The chevron 
shape is also seen in our other integrations, and indicates that eddies transport momentum 
equatorward (see Section 3.6). 

3.2. Effect of eccentricity 

Our models show that, at constant orbit-mean stellar flux, the dayside temperatures, 
day-night temperature differences, and wind speeds at periapse passage all increase with 
increasing orbital eccentricity. This is illustrated in Figure 7, which compares three 
simulations with an equilibrium temperature of 951 K ((F) = 185691 Wm~^), that vary 
in eccentricity from 0.0 (top row) to 0.75 (bottom row). Again, we plot each column 
on the same colorscale (zonal wind and temperature) for comparison. As eccentricity 
increases, peak temperatures increase from 1000 K to 1300 K. In addition, the variation 
in temperature from dayside to nightside increases; this strengthens the winds within the 
equatorial superrotating jet. Peak winds reach ^2500 ms~^ in the circular case, but increase 
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to ^5000 ms~^ in the high eccentricity case. 

As the eccentricity increases, the equatorial jet also narrows in latitude. In the case 
of highest eccentricity, the atmosphere develops westward midlatitude jets and eastward 
high-latitude jets. The narrowing of the equatorial jet arises from the fact that jet width is 
confined to the Rossby radius of deformation (Showman & Polvani 2011), the length scale 
over which the atmosphere adjusts to planetary-scale phenomena. The Rossby radius of 
deformation (Holton 2004) is defined at the equator as 



Ln=^-^. (6) 

where c is the gravity wave speed, and /? is the derivative of the Coriolis parameter, /, with 
respect to latitudinal distance — ^ — ^^^^^^^^ ^ where Qrot is the planetary rotation 
rate, is the latitude and Rp is the planetary radius. At the equator, /5 is simply 
Hence, at the equator Lji is 



V "^^rot 

Therefore, oc J — , and an increased rotation rate leads to a narrower deformation 
radius. From the Hut (1981) formulation, the combined effect of the larger orbital period 
and increase in eccentricity produces a shorter rotational period. This, in turn, increases the 
rotation rate, which reduces the deformation radius. For example, the ratio between Qrot 
for the e=0.0 case and the e=0.75 case is 6.3; hence, we would expect that the equatorial 
jet should be narrower by a factor of a/O ^ 2.5. This is indeed the case, as the e = 0.75 
case has an average jet width of ~ 40° latitude, while the e = 0.0 case has an average jet 
width of ~ 100° latitude. 



3.3. Effect of varying average stellar flux 



An increase in average stellar flux leads to increases in temperature and wind speeds, 
illustrated in Figure 8. Each case shown has an eccentricity of 0.25; these maps have 
independent colorscales to better show atmospheric structure. The planet that receives 
the highest average stellar flux throughout its orbit is hottest overall (bottom), with peak 
temperatures of 1300 K. Furthermore, because the planet with the highest (F) orbits at the 
smallest orbital distance, it has a faster rotation rate and hence a narrower equatorial jet; 
this comparison makes the effect of rotation on the jet structure much more apparent. 

Interestingly, the wind and temperature structure of the farthest case (a = 0.2019 AU, 
top row), is unlike the other three; the flow is fully westward at pressure levels above 1 mbar. 
Peak temperatures reach only 390 K, and no strong day-night temperature variations occur. 
This would suggest that a regime shift occurs between orbital distances of 0.08 and 0.20 
AU. The homogeneity in temperature is due to the planet's long radiative time constant 
and the slow rotation rate. The planet effectively transports heat both latitudinally and 
longitudinally, thus leading to small temperature variations. Exploring this regime shift in 
depth will be the subject of a future paper. 



3.4. Synchronous vs. pseudo-synchronous rotation 

Here we explore the sensitivity of the circulation to rotation rate by comparing 
models performed with the pseudo-synchronous rotation rate to otherwise identical models 
performed with the synchronous rotation rate. This is illustrated in Figure 9. Here, the 
rotation rate varies by almost a factor of 3. Both cases still have an equatorial superrotating 
jet, and the overall atmospheric structure is similar. However, there are subtle differences: 
in the synchronous case, the equatorial eastward winds and high-latitude westward winds 
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are higher in magnitude (right column) . The synchronous case also exhibits more latitudinal 
and longitudinal temperature variation, particularly on the dayside; this is more apparent 
at apoapse (bottom row). In the pseudo-synchronous case, the faster rotation rate serves to 
homogenize the temperature in longitude, presumably because the time for rotation to carry 
air parcels from day to night is shorter. The temperature difference from equator-to-pole is 
also larger, perhaps because the faster rotation rate inhibits equator-to-pole heat transport. 
Thus, while it is unlikely that an eccentric planet synchronously rotates its star, the 
distinction can lead to notable differences in the mean flow. 

3.5. Momentum budget and Equatorial superrotation 

In most of our simulations, the planet maintains an equatorial eastward jet throughout 
its orbit. Superrotation is common in three-dimensional circulation models of hot Jupiters 
(e.g. Showman & Guillot 2009, Showman et al. 2008, 2009, Menou and Rauscher 2009, 
Rauscher and Menou 2010, Cooper and Showman 2005,2006, Dobbs-Dixon & Lin 2008, 
2010, Heng et al. 2011), and is also seen within our own solar system on Venus, Titan, 
Jupiter and Saturn. While superrotation on solar system planets has been extensively 
studied, only recently have the mechanisms for superrotation on synchronously rotating 
exoplanets been identifled (Showman and Polvani 2011, hereafter SP2011). As discussed in 
their paper, the generation of superrotation requires up-gradient momentum transport that 
pumps angular momentum from outside the jet to within it. Only transport by waves and 
eddies can satisfy this requirement (Hide 1969). 

SP2011 demonstrate that for synchronous hot Jupiters on near-circular orbits such as 
HD 189733b and HD 209458b, equatorial superrotation is generated by standing, planetary 
scale Kelvin and Rossby waves at the equator and midlatitudes, respectively, that are 
themselves the dynamical response to the planet's large day-night heating contrast. Because 
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the Kelvin waves propagate eastward and the Rossby waves propogate westward, a phase 
tilt of the wind vectors emerges from northwest to southeast in the northern hemisphere, 
and southwest to northeast in the southern hemisphere, with an overall shape resembling a 
chevron pointing eastward. This pattern causes meridional eddy angular momentum fluxes 
that converge angular momentum onto the equator, generating equatorial superrotation. 

Despite the fact that not all our model integrations are synchronous, our results 
suggest that similar mechanisms are still at play. The chevron shape that SP2011 describe 
is exactly what we noted earlier in our plots of the wind and temperature profiles from 
Figure 5, suggesting that, as in SP2011, our models exhibit an equatorward flux of eddy 
angular momentum. This is demonstrated quantitatively in Figure 10, which shows the 
zonally averaged meridional flux of relative zonal eddy angular momentum, u'v' coscj)^ at 
snapshots throughout the orbit (as in Figures 5 and 6). Here, and represent the 
devation of the zonal and meridional winds, respectively, from their zonal averages. Positive 
values indicate a northward flux of eastward eddy angular momentum, and negative values 
indicate a southward flux of eastward eddy angular momentum. At apoapse, the small 
day-night forcing causes poleward transport of eddy angular momentum at latitudes greater 
than 50° but equatorward transport of eddy angular momentum equator of 50° latitude; 
the amplitudes are weak due to the low insolation at apoapse. As the planet moves closer 
to the star (middle row, left column), the increase in day-night forcing strengthens the 
equatorward flux of eastward eddy angular momentum at latitudes between 0° and 50° 
in both hemispheres. The momentum flux grows and persists until hours after periapse. 
Angular momentum thus converges onto the equator, and maintains the equatorial jet 
against westward accelerations caused by advection, friction, and Coriolis forces. 

These mechanisms also explain the difference in wind speeds between the synchronous 
and pseudo-synchronous cases in Figure 9. The synchronous case has a higher day-night 
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heating contrast, hence a higher equatorward flux of eddy momentum, leading to a stronger, 
narrower equatorial jet. Moreover, the equatorial jet in the pseudo-synchronous case 
extends to deeper pressures than that of the synchronous case; this is probably due to a 
change in momentum budget. 

4. Observational Implications 

Observations of transiting exoplanets in circular orbits are known to be shaped by the 
atmospheric circulation. As was flrst predicted by Showman & Guillot (2002) and later 
observationally confirmed by Knutson et al. (2007), an equatorial superrotating jet causes 
an eastward displacement of the hot spot from the substellar point under appropriate 
conditions. Hence, a lightcurve shows a peak in infrared flux minutes to hours before 
secondary eclipse. We expect such a situation for eccentric exoplanets as well. However, 
observations of eccentric planets can be complicated because the flux received by the planet 
varies throughout its orbit. Previous studies explored such situations, but have not used 
three-dimensional models (e.g. Langton and Laughlin 2008, Cowan and Agol 2011). We 
can expand on this study using our three-dimensional SPARC model results. 

We choose to generate lightcurves following the procedures described in Fortney 
et al. (2006) for each model integration at three values of the argument of periastron 
(cj), the angle between the radius vector to the ascending node and the periapse of the 
orbit. In particular, we choose orientations where transit occurs before periapse (at 
/ = — 90°;cl; = 180°), transit occurs at apoapse {uj = 270°), and transit occurs after 
periapse (at / = +90°; = 360°) (Figure 11). These values reflect the broad range of 
viewing geometries of transiting exoplanets: for example, HAT-P-2b (e = 0.52), HAT-P-17b 
(e = 0.35) and HD 97658b (e = 0.13) have co near 180°; WASP-8b (e = 0.31), HAT-P-31b 
(e = 0.25) and HD 80606b (e = 0.93) have u near 270°; and GJ 436b (e = 0.15), X0-3b 
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(e = 0.26), and HAT-P-llb (e = 0.20) have uo near 360° (see exoplanet.eu for a full list of 
values). Note that when transit occurs, the nightside is visible to Earth, while the dayside 
is visible around the time of secondary eclipse. 

4.1. Observational effects of varying uu 

Using the nominal case (e = 0.50, a = 0.0534 AU, T^q — 951 K), we illustrate the 
dependence of full-orbit lightcurves on the value of uj. Figure 12 shows three lightcurves 
of the same model integration, but at the differing values of uj described in the previous 
section. Hence, the differences in lightcurves results solely from the differing viewing 
geometries. The planet/star flux ratio is plotted as a function of Earth days from periapse 
in the warm Spitzer wavelength bands (3.6 and 4.5 /im), and JHK bands (1.26, 1.65, 
and 2.20 /am). The solid and dashed lines correspond to transit and secondary eclipse, 
respectively. For the lightcurve where uu — 180° (top panel), the peak in infrared (IR) flux 
occurs approximately 6 hours after periapse and before secondary eclipse. The timing of 
the peak is due to a combination of the radiative lag of the planet and the fact that the 
hottest regions rotate into an Earth-facing orientation hours after periapse passage. The 
dayside, which is most visible near secondary eclipse, is hottest after periapse. Additionally, 
temperatures reach peak values at times close to when those hottest regions are facing 
Earth. This helps explain the high amplitude of the lightcurve peak in the uj = 180° case 
relative to the other two cases. These effects can be visualized by plotting the temperature 
near the photosphere along the equator in a polar projection (Figure 13). Snapshots 
throughout the orbit plot the temperature in a polar projection as a function of pressure 
and orbital position at a latitudinal slice near the equator; the solid and dashed lines denote 
the substellar longitude and earth-facing longitude, respectively. Based on the temperatures 
at the earth-facing longitude, it is apparent that the peak in IR flux should occur after 
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periapse and before secondary eclipse (see / = 12° and 72°); this is indeed the case. 

For the same model but instead with uo = 270°, the peak in IR flux occurs 1 hour after 
eclipse and periapse. In this case, the dayside is visible near periapse, and therefore we 
are seeing the planet increase in temperature (and IR flux) during its closest approach to 
the star. However, the hottest regions are displaced eastward of the substellar point, and 
therefore the hottest regions face Earth signiflcantly before periapse passage (Figure 13 at 
/ 72°). Thus, by the time the peak temperatures are reached near and after periapse, 
the hottest regions are already rotating out of view from Earth. This explains the lower 
amplitude of the lightcurve peak as compared to the peak of the uo = 180° case. The IR 
peak decreases as the planet moves further away from the star along our line of sight. A 
second ramp-up in flux occurs approximately 1 day after periapse/eclipse; this is due to the 
hottest regions rotating back into view as the planet moves toward apoapse. 

In the case where uj — 360°, two distinct peaks are present - a smaller peak that occurs 
after secondary eclipse and ^3-4 hours before periapse, and a second higher peak that 
occurs after transit. The flrst peak occurs as some of the dayside is visible as the planet is 
heated approaching periapse. The peak decreases as the dayside rotates out of view. For 
this geometry, the hottest regions are facing away from Earth when they reach their peak 
temperatures (Figure 13). The second peak occurs after transit as the hottest regions again 
rotate into view. At this point in the orbit, peak temperatures have already occurred, and 
therefore this peak is lower in amplitude than the other two lightcurves. However, because 
the planet is hotter than it was after secondary eclipse, this second peak has a higher 
amplitude than the flrst one. 



4.2. High-e model integrations at varying (F) 



If an eccentric planet has a large enough day-night temperature difference and a 
sufficiently rapid rotation rate, its resultant lightcurve can show a periodic rise and fall in flux 
throughout its orbit. The lightcurves shown in Figure 14 are two e = 0.75 model integrations 
that differ in (F); the top panel shows a case where (F) = 468183 W m"^ {T^q = 1199 
K) and a = 0.0385 AU, while the bottom panel is a case where (F) = 185691 W m~^ 
{Teq = 951 K)and a = 0.0611 AU. For both cases, lightcurves are shown for the uj = 180° 
geometry (see Figure 11). The differences in (F) lead to differences in lightcurve amplitude 
in the two models; the closer, hotter planet exhibits a higher flux amplitude. Not 
surprisingly the timing of peak IR flux for both models resembles the ou — 180° case in 
Figure 12. The peak in each case occurs after periapse and before eclipse, due to the 
combined effects of the thermal lag of the planet and geometry of the system. Additionally, 
both lightcurves exhibit a quasi-periodic rise and fall of IR flux days after periapse passage. 
This phenomenon, called "ringing" , was also seen in models by Langton & Laughlin (2008) 
and Cowan & Agol (2011), caused by the hottest point of the planet rotating in and out of 
view from Earth. For ringing to occur, the temperature difference from dayside to nightside 
must be large and it must survive for multiple planetary rotation periods. In the models 
shown in Figure 14, this is aided by the rapid pseudo-synchronous rotation rate associated 
with an eccentricity of 0.75 (Table 2). For the top case, the day-night temperature difference 
is over ^ 1300 K, while the bottom case varies by ^ 750 K (Figure 7). The lightcurves 
from Figure 12 do not exhibit ringing because the day-night temperature difference is low 
(- 300 -400 K). 

Because the radiative and dynamical timescales in the models by Cowan & Agol (2011) 
are calculated by a simple scaling with temperature, the ringing seen in their simulated 
lightcurves occurs with a period equal to the planet's assumed solid-body rotation period 
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in the inertial frame of the winds. In reahty, an eccentric exoplanet approaching apoapse 
would become increasingly cold, the winds would weaken and decrease, and hence the 
ringing should be non-periodic. Indeed, for our three dimensional models, the ringing in 
each case has varying periods - the top case has a ringing period ranging from 0.3 to 0.4 
days (Trot — 0.35 days) while the bottom panel has a period ranging 0.7-0.8 days (Trot — 0.7 
days) . 

Despite these differences, this study and the studies by Langton & Laughlin (2008) and 
Cowan & Agol (2011) illustrate the influence of atmospheric dynamics on the observations 
of eccentric exoplanets. The lightcurve features described here could be seen in future 
follow-up observations of eccentric exoplanets, and observers can use the conclusions drawn 
above to relate the observations to the planet's atmospheric structure. For example, if 
ringing is present in an IR lightcurve, the planet must have a strong day-night temperature 
assymmetry, and hence we would expect the planet to have strong superrotation at the 
equator. This is probably the case for highly eccentric (e > 0.5), short-period exoplanets 
such as HD 17156b and HD 80606b. If ringing were observed, it would also allow an 
estimate of the sum of the rotation speed and the superrotation speed, thereby placing 
constraints on both wind speeds and rotation rates of eccentric hot Jupiters. The timing of 
peak IR flux relative to transit and secondary eclipse, and also relative to orbital position, 
can help to constrain the circulation further. 



5. Conclusions 

We present three-dimensional circulation models coupled with a two-stream, non-grey 
radiative transfer scheme for a number of theoretical eccentric hot Jupiters. We have shown 
that as in published models with zero eccentricity, our high-eccentricity circulation models 
are dominated by eastward flow at photospheric levels which cause an eastward displacement 



-20- 



of the hottest regions from the substellar point. The rapid rotation rates associated with 
pseudo-synchronization at high eccentricity lead to a small Rossby deformation radius and 
in some cases multiple jets in the atmosphere. Global-mean temperatures and day-night 
temperature differences peak not at periapse but several hours afterward. 

Furthermore, we show that equatorial superrotation is generated and maintained by 
eddies formed by the strong day-night heating contrast, which induce a flux of momentum 
from midlatitudes to the equator. The eddy magnitudes and momentum fluxes peak just 
after periapse passage leading to variations in the zonal-mean flow throughout the orbit. 

Lastly, we have shown that the spatial and temporal variations of the wind and 
temperature structure, as well as the orbital viewing geometry of the system with respect 
to Earth, can affect the time and amplitude of peak IR flux seen in full-orbit lightcurves. 
Depending on the viewing geometry of the orbit relative to Earth, we flnd that peaks in IR 
flux that either lead or lag periapse are possible; in all cases, a combination of temporal 
effects (temperatures changing over time) and geometric effects (hot spots rotating into or 
out of view) are important in controlling the timing and amplitude of the flux peaks. In 
cases where the day-night temperature contrast is large and the rotational period is short, 
the lightcurve can also exhibit "ringing" in flux as the hottest region of the planet rotates 
in and out of view. This ringing is non-periodic, due to the variation in stellar heating as a 
function of distance. 

This work was supported by NASA Origins and Planetary Atmospheres grants to 
APS. T.K. also acknowledges support from the Harriet P. Jenkins Pre-Doctoral Fellowship 
Program (JPFP). 



- 21 - 



A. Validation tests for updated radiative transfer scheme 

To validate our updated radiative transfer scheme using 11 frequency bins, we 
conducted the following tests: 

1. First, we tested the sensitivity of the number of frequency bins in a one- 
dimensional radiative-convective model of HD 189733b (a = 0.0313 AU, e = 0.0, 
Rp — 8.2396 X lO'' m, g = 21.4 ms~^). Figure 15 compares the globally averaged 
pressure-temperature profiles of models run with 11 frequency bins, 30 frequency bins, 
and 196 frequency bins at (from left to right) Ix, 2x, 4x, and 8x the average stellar 
flux of HD 189733b. For each group of models, the 11-bin proflle agrees well with the 
30- and 196-bin cases from 1 /ibar to 1 bar; the temperature varies by only tens of K. 
Below 10 bars, the 30- and 11-bin profiles differ from the 196-bin cases by up to 50 K, 
due to the loss of resolution at low wavenumbers. In the case of the highest average 
flux, the large amount of heating leads to a temperature inversion that is not fully 
captured by the 11-bin model, but can be improved with a better initial guess at the 
P-T profile. 

2. After testing the new RT scheme for one-dimensional models, we ran full three- 
dimensional simulations of HD 189733b at its nominal average stellar fiux using 
the SPARC/MITgcm. Figure 16 shows the globally-averaged PT profiles (top row), 
zonal-mean zonal wind (middle row) and wind/temperature profiles at 30 mbar 
(bottom row) for the 30-bin (left column) and 11-bin (right column) RT setups. The 
globally-averaged profiles are nearly identical, with any temperature differences within 
a few K. The zonal wind plots also show good agreement, with similar peak speeds 
(> 2600 ms~^), equatorial jet width (^ 60°), and jet level (< 100 bar). At 30 mbar, 
the 30- and 11-bin models share similar dayside-nightside temperature structure, with 
a hot spot (1100 K) eastward of the substellar point, and two colder regions on the 
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nightside in the mid-latitudes. 

Both tests show that the bulk circulation and temperature structure is retained in the 
transition from 30 to 11 frequency bins; hence, we proceed to use the SPARC/MITgcm 
with the new scheme for these and future model simulations. 
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Table 1. Bounding wavelengths of the frequency bins used in the radiative transfer 

calculations. 



Wavelength (/xm) 


Wavelength (/xm) 


324.68 


20.00 


20.00 


8.70 


8.70 


4.40 


4.40 


3.50 


3.50 


2.50 


2.50 


2.02 


2.02 


1.32 


1.32 


0.85 


0.85 


0.61 


0.61 


0.42 


0.42 


0.26 



Table 2. List of model integrations. 



. (W m-2) 


Tec (K) 


a (AU) 


e 


Tp (AU) 


Ta (AU) 


Torb (s) 


Trot (s) 


Torb /Trot 


Q (s- 






468183 


1199 


0.0313 


0.00 


0.0313 


0.0313 


1.917 X 


10^ 


1.917 X 


10^ 


1.0 


3.27 X 


10- 


-5 


468183 


1199 


0.0318 


0.25 


0.0239 


0.0398 


1.963 X 


10^ 


1.421 X 


10^ 


1.4 


4.42 X 


10- 


-5 


468183 


1199 


0.0385 


0.75 


0.00963 


0.0674 


2.615 X 


10^ 


3.010 X 


10"^ 


8.7 


2.09 X 


10- 


-4 


185691 


951 


0.0497 


0.00 


0.0497 


0.0497 


3.834 X 


10^ 


3.834 X 


10^ 


1.0 


1.64 X 


10- 


-5 


185691 


951 


0.0505 


0.25 


0.0379 


0.0631 


3.927 X 


10^ 


2.843 X 


10^ 


1.4 


2.21 X 


10- 


-5 


185691 


951 


0.0534 


0.50 


0.0267 


0.0801 


4.270 X 


10^ 


1.522 X 


10^ 


2.8 


4.13 X 


10- 


-5 


185691 


951 


0.0611 


0.75 


0.0153 


0.1069 


5.226 X 


10^ 


6.016 X 


10"^ 


8.7 


1.04 X 


10- 


-4 


73680 


755 


0.0789 


0.00 


0.0789 


0.0789 


7.668 X 


10^ 


7.668 X 


10^ 


1.0 


8.19 X 


10- 


-6 


73680 


755 


0.0802 


0.25 


0.0602 


0.1003 


7.858 X 


10^ 


5.689 X 


10^ 


1.4 


1.10 X 


10- 


-5 


73680 


755 


0.0848 


0.50 


0.0424 


0.1272 


8.544 X 


10^ 


3.046 X 


10^ 


2.8 


2.06 X 


10- 


-5 


11617 


476 


0.1987 


0.00 


0.1987 


0.1987 


3.067 X 


10^ 


3.067 X 


10^ 


1.0 


2.05 X 


10- 


-5 


11617 


476 


0.2019 


0.25 


0.1514 


0.2524 


3.141 X 


10^ 


2.274 X 


10^ 


1.4 


2.76 X 


10- 


-5 


11617 


476 


0.2135 


0.50 


0.1068 


0.3202 


3.416 X 


10^ 


1.218 X 


10^ 


2.8 


5.16 X 


10- 


-6 


11617 


476 


0.2443 


0.75 


0.0611 


0.4275 


4.181 X 


10^ 


4.813 X 


10^ 


8.7 


1.31 X 


10- 


-5 
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detected exoplanet 
468183 Wm"^ 
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Fig. 1. — Model integrations from this study plotted against a portion of the current exo- 
planet population (dots). The green right-pointing triangles, red down-pointing triangles, 
blue left-pointing triangles and magenta up-pointing triangles correspond to average stellar 
fluxes of 468183, 185691, 73680 and 11617 Wm-^, respectively. See Table 2 and text for 
more details. 
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Temperature (K) 

Fig. 2. — Initial 1-D pressure-temperature profiles assigned to each column of the atmosphere 
at the start of each simulation, at the four different values of average stellar flux. These 
profiles were calculated assuming radiative equilibrium, and use the same opacities as the 
full 3D model. 
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e= 0.50, a= 0.0534 AU 
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Fig. 3. — Root-mean squared velocity plotted as function of pressure and simulation time 
for one of our model integrations. The plot shows a case where Teq = 951 K, e = 0.50, and 
a = 0.0534 AU. RMS velocity was calculated at each pressure level from wind speeds output 
every 100 days. 
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Fig. 4. — Dynamical response of a model with T^q = 951 K, e = 0.50, and a = 0.0534 AU. 
Top: Global, horizontally averaged temperature versus pressure and time throughout the full 
orbit; zero is defined as periapse passage. The upper axis gives the true anomaly throughout 
the orbit. Bottom: zonally- averaged zonal (east- west) winds versus pressure and latitude for 
the same model, time- averaged over a full orbit. 
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Fig. 5. — Snapshots of the wind and temperature structure at 30 mbar of a planet with 
T^q = 951 K, e = 0.50, and a = 0.0534 AU. To trace the behavior throughout one complete 
orbit, the different panels should be foUoweed counterclockwise from the the top panel. 
Snapshots are taken from apoapse (true anomaly / = —180.0°, top panel) to hours after 
periapse (/ = 119° , follow arrows) and back again to apoapse. Each panel is plotted on the 
same temperature colorscale. Note that peak temperatures lag periapse passage. Periapse 
corresponds to / = 0° and apoapse to 180°. The vertical line corresponds to substellar 
longitude. 





Fig. 6. — Snapshots of the zonal-mean zonal wind for the same integration described in 
Figure 5 {T^q = 951 K, e = 0.50, and a = 0.0534 AU). Snapshots are again taken from 
apoapse (/ = —180.0°) to hours after periapse (/ = 119°) and back again to apoapse. Each 
panel is plotted on the same temperature colorscale. Note that the equatorial superrotating 
jet is persistent throughout the orbit, with peak speeds exceeding 5000 ms~^. 
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e= 0.00 




e= 0.50 




Fig. 7. — Comparison of simulations with T^q = 951 K with increasing eccentricity, from 
e = 0.0 (top row) to e = 0.75 (bottom row). Shown are plots of the zonal-mean zonal wind 
averaged over a full orbit (left column) and snapshots of the wind and temperature profiles 
at 30 mbar (right column) . The snapshots of the wind and temperature profiles in the right 
column were taken at periapse passage. The vertical bars in the right column denote the 
longitude of the substellar point. 



- 35 - 



a= 0.2019 AU 




-80 -60 -40 -20 20 40 60 80 (m s"^) 
latitude (°) 



a= 0.0802 AU 




-80 -60 -40 -20 20 40 60 80 (m s"^) 
latitude (°) 



a= 0.0505 AU 




-80 -60 -40 -20 20 40 60 80 (m s"^) 
latitude (°) 



3= 0.031 8 AU 




-80 -60 -40 -20 20 40 60 80 (m 3"^ 
latitude (°) 



a= 0.2019 AU 




-150 -100 -50 50 100 150 (K) 
longitude (°) 



a= 0.0802 AU 




-150 -100 -50 50 100 150 (K) 
lonqitude (°) 



a= 0.0505 AU 




-150 -100 -50 50 100 150 (K) 
lonqitude (°) 



a= 0.0318 AU 




-150 -100 -50 50 100 150 (K) 
longitude (°) 



Fig. 8. — Comparison of simulations with e = 0.25 with decreasing average steUar flux, 
from bottom to top panels. Shown are plots of the orbit- averaged zonal-mean zonal wind 
(left column) and the wind and temperature proflles at 30 mbar at periapse (right column). 
Each plot has an independent colorscale. The vertical bars in the right column denote the 
substellar longitude. 



- 36 - 



e = 0.50, pseudo-synchronous rotation 



synchronous rotation 




-150 -100 -50 50 100 150 
longitude (°) 
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Fig. 9. — Comparison of pseudo-synchronous (Trot = 3.046 x 10^ s) versus synchronous 
(Trot = 8.544 X 10^ s) rotation, for a case with T^q = 755 K, e = 0.50, and a = 0.0848 AU. 
Wind and temperature snapshots are taken at periapse (middle row) and apoapse (bottom 
row). The vertical bar denotes substellar longitude. 
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Fig. 10. — Eddy momentum flux, u'v' coscj)^ throughout the orbit of a model where T^q = 951 
K, e = 0.5, and a = 0.0534 AU. Positive values indicate a northward flux of eastward eddy 
angular momentum, while negative values indicate a southward flux of eastward eddy angular 
momentum. The maximum momentum flux occurs hours after periapse passage (/ = 72°). 
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Fig. 11. — Differing orbital viewing geometries explored in this study: u = 180°, co = 270°, 
cj = 360°. 



- 39 - 



CO 



c 

CO 



0.0015 


. cj = 180 deg 




' ' ' 1 1 

4.5/i.m 
3.6/i.m 
K 

\ ^ 


0.0010 






\ ^ 


0.0005 








0.0000 









-2-1012 
Time after periapse passage (days) 



0.0015 



OS 
DC 

§ 0.0010 



OS 

I 0.0005 

Q_ 



0.0000 





1 1 1 1 1 1 






4.5/Lim 




^ cj = 270 deg . 


3.6/im 






\ ^ 






\ ^ 























-2-1012 
Time after periapse passage (days) 



0.0015 



DC 

§ 0.001 Oi 



I 0.0005 [ 

CL 



0.0000 













4.5/zm 


cj = 560 deg 




3.6yLim 






K 






H 






J 












-2 -1 





1 2 



Time after periapse passage (days) 



Fig. 12. — Three lightcurves of the IR flux (expressed as a planet/star flux ratio) versus time 
throughout the fuU orbit for a single model integration where T^q = 951 K, e = 0.50, and 
a = 0.0534 AU. The flux is plotted in warm Spitzer bandpasses, 3.6 and 4.5 /xm, and JHK 
bands (1.26, 1.65, 2.20 /xm, respectively). These three cases only differ by their value of cj, 
shown in Figure 11. The solid lines denote transit, while the dotted lines denote secondary 
eclipse. 



^40- 




Fig. 13. — Plots of temperature versus longitude and pressure in equatorial slices throughout 
one full orbit, for our nominal model with T^q = 951 K, e = 0.50, and a = 0.0534 AU. Each 
annulus shows temperature versus longitude and log-pressure for a snapshot at a particular 
orbital position (denoted by the value of /) as viewed from looking down from over the 
planet's north pole. Pressures between 30 and 43 mbar are displayed in each annulus. 
The different annuli mark the time evolution throughout one complete orbit, starting from 
apoapse for the top annulus, proceeding to peripase, and then back to apoapse (follow arrows 
moving counterclockwise around the plot). The solid line indicates the longitude of the 
substellar point while the dashed line denotes the Earth-facing longitude for the particular 
case oi u) = 180°. For the cases of uj = 270 and 360 degrees, exactly the same plot is valid, 
except that the sub-Earth longitude would be aiming upward for uj = 270° and right for 
uu = 360°. 
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Fig. 14. — Lightcurves for two different e = 0.75 model integrations: Tgg = 1199 K, a = 
0.0385 AU (top), and Teq = 951 K, a = 0.0611 AU (bottom), both with co = 180°. The sohd 
Hne denotes transit, while the dotted line denotes secondary eclipse. 
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Fig. 15. — Globally-averaged pressure temperature profiles from 1-D radiative-convective 
radiative transfer models of HD 189733b with 11 frequency bins (red profiles), 30 frequency 
bins (green profiles) and 196 frequency bins (black profiles). From left to right, the profiles 
vary in stellar fiux from the nominal value, as well as 2x, 4x and 8x higher fiux. The 8x 
higher fiux case has a temperature inversion. 
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Fig. 16. — Globally-averaged pressure-temperature profiles (top row), zonal-mean zonal wind 
(middle row) and the wind/temperature profiles at 30 mbar (bottom row) for the 30-bin (left 
column) and 11-bin (right column) model integrations of HD 189733b. 



